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Use  of  an  ejector  is  an  effective  way  to  increase  the  thrust 
produced  by  a  jet.  In  this  thesis  project  an  axisymmetric 
ejector  concept  which  had  been  previously  explored  by 
experiment  was  numerically  modeled.  An  existing 
axisymmetric,  internal  flow  code  based  on  the  explicit 
MacCormack  method  was  modified  to  incorporate  primary  nozzle 
structure  and  flow  injection  within  the  flowtield.  Results 
were  compared  qualitatively  and  quantitatively  with 
experimental  results  to  verify  the  validity  of  the  model. 
Internal  flow  structure,  difficult  to  obtain  in  experiment, 
is  easily  examined.  This  code  may  be  used  for  parametric 
analysis  of  such  ejector  performance  parameters  as  primary 
nozzle  location,  flow  injection  angle,  diffuser  area  ratio, 
and  inlet  geometry  to  optimize  future  hardware 
configurations.  , 
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NUMERICAL  ANALYSIS  OF  AN  AXISYMMETRIC  THRUST 


AUGMENTING  EJECTOR 


I  Introduction 


Background 

Ejectors  provide  an  attractive  solution  to  a  problem 
which  has  plagued  designers  of  jet  powered  Vertical/Short 
Takeoff  and  Landing  (V/STOL)  aircraft  since  they  were 
conceived.  The  problem  being,  that  in  order  to  get  thrust 
greater  than  the  weight  of  the  vehicle  such  that  vertical 
takeoff  or  landing  can  be  achieved,  an  aircraft  must  be 
equipped  with  a  powerplant  which  is  larger  (sometimes 
substantially)  than  that  which  is  needed  for  its  primary 
mission.  Ejectors,  due  to  their  thrust  augmenting 
characteristics,  provide  a  means  for  achieving  thrust  to 
weight  ratios  greater  than  unity  with  a  more  reasonably 
sized  powerplant. 

In  an  ejector,  a  jet  flow  is  directed  into  an  open  duct 
of  somewhat  larger  cross-sectional  area.  A  secondary  flow 
is  induced  through  the  duct  by  viscous  entrainment  (Figure 
1).  A  simple  analysis  de\ sloped  by  Von  Karman  (18:461-62) 
shows  that  the  thrust  theoretically  achievable  by  such  a 
process  is  up  to  twice  that  of  the  primary  nozzle  alone 
(Figure  2). 
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The  thrust  augmentation  ratio  (0)  is  defined  as  the 
ratio  of  the  total  thrust  of  the  ejector  (F)  to  the 
isentropic  thrust  of  the  primary  nozzle  alone  (assuming 
incompressible  flow): 

0  =  1  /  F rtoz  isen 

<p  =  Q  v\  A  /  Q  vj  Aj 

<P  =  (Ve/V,)2  /  a  (1) 

where  cx  is  the  ratio  of  the  area  of  the  jet  (Aj)  to  the  area 
of  the  mixing  chamber  (A) 

a  =  Aj/A 

and  Vh/Vj  is  the  ratio  of  the  flow  velocity  following 
complete  mixing  (Vb)  to  the  flow  velocity  of  the  primary 
nozzle  (Vj)  found  by  applying  continuity  f’-om  the  jet  exit 
plane  to  the  exit  of  the  ejector: 


V\  (  A  -  Aj  )  +  ^  Aj  =  VE  A 


(2) 


and  applying  conservation  of  momentum: 

A  (P\  -  Pn)  =  g  A  V\-q  Aj  Vf  -  q  (A -A,)  V\  (  3  ) 

Applying  Bernoulli's  equation  from  ahead  of  the  inlet  to  the 
jet  exit  plane: 


/J„  =  Px  +  /z  Q  V\ 


(4) 


gives  three  equations  in  three  unknowns,  Vb,  Vi,  and  Pi. 
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Eliminating  Vi  and  Pi,  the  velocity  ratio,  Ve/Vj ,  is  found 
to  be : 


I  ~a  (  1  -  2a  )  +  (  2a  -  6a 2  +  6a''  -  2a4  ] 

(  !  -  2a  +  2a2  ) 


(5) 

(7:28) 


The  thrust  augmentation  achieved  may  be  further 
increased  through  the  use  of  a  diffuser  (Figure  3). 


Through  a  similar  analysis,  the  modified  equation  for  a 
diffuser  equipped  ejector  is  found  to  be: 


2afi  ( 1  -  2a) 

2a  -  3a2 

1  -  2a  +  a2(  1  +  ft2) 

1  -  2a  +  a2(\  +  ft2) 

0  (6) 
(9:282) 


([/,/!/)-  +  (Vy/Vj) 


with  0  now  given  by  (assuming  incompressible  flow): 


<p  =  Q  V\  A,  /  q  V]  A, 

<P  =  (Vi-JV,)1  (A,./ A)  [A /A,) 


<P  =  (VE/V>r  05/a) 


(7) 


where : 


fi  =  At! A 

with  Ve  and  Ae  as  defined  in  Figure  3  (9:282-3). 

Plotting  the  thrust  augmentation  ratio  shows  that  for 
every  primary  nozzle  to  mixing  chamber  area  ratio  (Aj/A) 
there  is  an  optimum  diffuser  area  ratio  (Ae/A)  which 
provides  the  maximum  thrust  augmentation  (Figure  4). 

Within  the  practical  limit  of  area  ratios,  it  can  be  seen 
that  thrust  augmentations  of  up  to  3  times  the  primary  jet 
thrust  are  theoretically  achievable. 

These  theoretical  values,  however,  assume  that  the 
ejector  is  long  enough  that  complete  mixing  of  the  primary 
and  secondary  flows  takes  place.  In  practice,  the  size  of 
the  ejector  is  limited  by  aircraft  installation  constraints 
and  frictional  effects,  which  lead  to  separations  in  long 
diffusers.  In  such  short  ejectors  mixing  is  incomplete  and 
coanda  nozzles,  in  which  the  jet  flow  is  injected  along  the 
wall  (Figure  5),  have  proven  to  be  more  effective  due  to  the 
greater  mixing  surface  area  of  the  primary  jet.  (9:286) 
Coanda  nozzles  provide  an  additional  benefit  for 
ejectors  equipped  with  a  diffuser  in  that  the  high  energy 
flow  along  the  wall  helps  to  prevent  separation  in  the 
diffuser  allowing  greater  area  ratios  to  be  achieved. 
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Fig.  4.  Thrust  Augmentation  as  a  Function 

of  cx-  and  J3  .  (9:284) 


Fig.  5.  Ejector  with  Coanda  Nozzles  (9:285) 
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In  a  previous  study,  Kedem,  (6)  experimentally  inves¬ 
tigated  a  two  dimensional  ejector  with  center  and  coanda 
nozzles  and  found  that  coanda  nozzles  did  provide  better 
performance  (6:51).  Reznick  (11)  and  later  Unnever  (17) 
expanded  on  these  results  by  testing  both  two  dimensional 
and  circular  ejectors  with  coanda  nozzles.  The  circular 
ejectors  provided  the  best  performance,  achieving  a  thrust 
augmentation  ratio  of  2.09.  This  was  accomplished  through 
the  use  of  discrete  jets  located  about  the  periphery  of  the 
inlet  which  are  referred  to  as  'hook  nozzles'  because  of 
their  hook-like  shape  (11:44).  This  approach  was  further 
investigated  by  Lewis  (8)  and  later  Uhuad  (16).  Hardware 
limitations  necessarily  constrained  the  various 
configurations  which  could  be  tested.  It  is  desirable  to 
have  a  computational  model  of  a  circular  ejector  which  is 
free  of  such  limitations.  With  this  model  a  full  range  of 
inlet  and  diffuser  geometries,  nozzle  sizes  and  positions, 
and  nozzle  plenum  pressures  can  be  investigated.  The 
computational  solutions  also  provide  complete  data  on  the 
behavior  of  the  flow  within  the  ejector,  which  is  difficult, 
if  not  impossible,  to  obtain  experimentally.  This  model 
has,  however,  been  constrained  to  an  axisymmetric  case  to 
reduce  complexity  and  required  computational  power.  This 
eliminates  the  possibility  of  modeling  the  hook  nozzles  used 
to  achieve  the  high  thrust  augmentation  observed  by  Reznick. 
This  model  is  therefore  constrained  to  the  case  of  an 
annular  nozzle,  which  was  also  one  of  the  configurations 
tested  by  both  Reznick  and  Unnever. 
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Purpose  and  Objectives 


The  primary  objectives  of  this  study  are: 

1.  Develop  a  numerical  algorithm  that  models  a  coanda 
flow  axisymmetric  ejector  of  the  type  tested  by  Reznick  and 
Unnever . 

2.  Validate  the  model  by  demonstrating  agreement  with 
experimental  results. 

3.  Provide  data  visualization  of  internal  flow 
characteristics . 

4.  Provide  sufficient  flexibility  to  accomplish 
secondary  objectives. 

5.  Begin  preliminary  investigation  of  secondary 
obj  ect ives . 

The  secondary  objectives  are: 

1.  Investigate  the  geometrical  and  fluid  injection 
parameters  which  effect  ejector  performance.  The  parameters 
which  may  be  investigated  include: 

a)  Fluid  injection  angle 

b)  Nozzle  position 

c)  Nozzle  area  ratio 

d)  Diffuser  area  ratio 

e)  Diffuser  geometry 

f)  Pressure  ratio  (Po  no*  /  P«) 

2.  From  this  parametric  analysis,  define  new  ejector 
configurations  which  may  merit  further  experimental 
investigation . 
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Scope 


A  generic  axisymmetr ic ,  internal,  Navier-Stokes  flow 
solver  based  on  the  explicit  MacCormack  method  provided  by 
WRDC/FIMM  was  modified  to  incorporate  primary  nozzle 
structure  and  flow  injection  within  the  flowfield  (program 
JSIAXK,  App  B).  A  diffuser  geometry,  nozzle  pressure  ratio, 
and  nozzle  exit  area  and  location  were  chosen  for  each  run. 
The  program  was  then  run  to  a  steady-state  solution  for  a 
given  flow  injection  angle  (oli).  The  injection  angle  was 
then  varied  and  the  run  continued  to  new  steady-state 
solutions  to  collect  data  oncti  -  variation  of  augmentation. 
Additional  runs  were  made  with  different  nozzle  locations 
holding  pressure  ratio,  nozzle  area,  and  nozzle  distance 
from  the  wall  constant  to  provide  data  on  nozzle  location 
effects  on  augmentation.  Also,  runs  were  made  with 
different  nozzle  areas  for  the  same  location  to  provide  data 
on  area  ratio  effects.  Runs  with  different  geometries  and 
pressure  ratios  could  provide  data  on  variations  with  these 
parameters . 

The  achieved  thrust  was  calculated  by  a  control  volume 
analysis  (Figure  6).  This  analysis  results  in  the  equation: 

F  =  j ^  QU2dA  (8) 

where  the  right  side  of  the  equation  is  the  integral  of  the 
x  -  momentum  at  the  exit  (5:12-13).  The  achieved  thrust  was 
compared  to  the  isentropic  thrust  of  the  primary  nozzle 
alone  to  calculate  the  thrust  augmentation  ratio  as  in 
equation  (5).  Thrust  augmentation  ratios  were  compared 
directly  with  experimental  data  for  similar  configurations 
and  the  differences  analyzed.  Exit  plane  velocity 
distributions  and  variation  of  thrust  augmentation  with  ot-i 
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and  position  were  compared  qualitatively  with  available 
experimental  results. 


Fig.  6.  Control  Volume 
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The  geometry  of  the  axisymmetric  ejector  which  was 
investigated  by  Reznick  and  Unnever  and  was  modeled  in  this 
study  is  portrayed  in  Figure  7.  The  mixing  chamber  was  4.4 
inches  in  diameter  and  3  inches  in  length,  the  inlet  lip  had 
a  2-inch  radius  of  curvature,  and  the  multi-section  diffuser 
could  be  changed  to  alter  the  diffuser  angles  (V i.'Vz,  ...) 
and  lengths  (Li,  L2,  ...),  up  to  three  sections,  and  thus 

the  area  ratio  Ae/A.  The  annular  exit  of  the  nozzle  had  a 
gap  of  0.065  inches  and  a  radius  of  3.0025  inches  and  could 
be  positioned  at  any  distance  (Li)  from  the  ejector  throat 
(constant  area  mixing  chamber).  (17:52) 


Fig.  7.  Ejector  Geometry 
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The  computational  grid  was  generated  us 
EJECGRD  (Appendix  A-l)  provided  by  WRDC/FIMM 
for  the  required  geometry.  Taking  advantage 
e  grid  models  a  single  radial/ax 
ne  to  the  wall  (Figure  8). 


Plane  of  Modeled  Flow 


ial)  grid  lines  are  defined  as 


circular  arcs  perpendicular  to  the 


all  and  ce 


For  the  mixing  area  and  diffuser,  by  parabolas 
perpendicular  to  the  wall  and  centerline  (which 
degenerate  to  straight  lines  for  the  mixing  area). 

Spacing  between  horizontal  grid  lines  was  varied 
exponentially  to  refine  the  grid  near  the  wall  where  large 
gradients  are  found  due  to  the  jet.  The  grid  was  also 
refined  near  the  centerline  to  minimize  the  effects  of  the 
boundary  conditions.  Spacing  between  vertical  grid  lines 
was  varied  exponentially  to  refine  the  grid  in  the  vicinity 
of  the  inlet  (and  the  jet  exit)  to  help  offset  the  tendency 
of  the  griu  to  spread  out  due  to  the  increasing  radii  of  the 
circular  arcs  in  the  inlet  and  also  to  concentrate  on  the 
larger  gradients  upstream.  The  exponential  spacing  was 
accomplished  using  the  program  DEHNEN  (Appendix  A-2)  also 
provided  by  WRDC/FIMM  which  was  used  to  produce  vertical  and 
horizontal  grid  spacing  tables,  EJECVRT  and  EJECHRZ 
(Appendices  A-3,~4).  Once  a  nozzle  location  was  identified, 
the  grid  was  refined  both  vertically  and  horizontally  in  the 
area  of  the  nozzle  exit.  An  example  of  the  grid  is  shown  in 
Figure  9.  Note  that  only  every  fourth  grid  line  has  been 
drawn  for  the  purpose  of  clarity. 

Flowfield  Initialization 

The  program  EJECGRD  also  calculates  an  initial 
flowfield  based  on  a  simple  one-dimensional  flow  area  ratio 
formula  given  a  value  of  the  critical  area  (A*),  the  duct 
area  through  which  the  flow  would  theoretically  be  sonic. 
This  critical  area  was  chosen  to  provide  a  reference  Mach 
number  at  the  diffuser  exit  of  0.1  (111.7  ft/sec).  This 
program  was  modified  to  overlay  a  high-speed  jet  flow  over 
the  low  speed  duct  flow  calculated  by  this  procedure.  This 
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was  necessary  because  the  formula  used  assumes  a  mass  flow 
balance  between  the  beginning  and  end  of  the  duct.  If  the 
jet  flow  is  not  accounted  for,  this  provides  an  extremely 
pior  initial  flowfield  which  can  lead  to  extended  processing 
time  or  even  cause  the  solution  to  'blow  up',  or  diverge,  as 
the  flowfield  tries  to  adapt  to  the  increased  mass  flow 
injected  by  the  primary  nozzle. 


F  Iqh  _S.fl.lv  er 

The  flow  solving  program  JSIAXK  (Appendix  B)  solves  the 
conservative  form  of  the  Navier-Stokes  equations  using  the 
explicit  MacCormack  scheme.  The  conservative  Navier-Stokes 
equations  in  axisymmetric  coordinates  are  (3:3): 

Equation  of  Mass 

5te,+  »W">  -0  (9) 

Equation  of  Axial  Motion 

+  '--(0"2)  +  U/r)  y( Quvr )  =  +  (1/r)  (rrxr)  /  ly  » 

<>t  dx  dr  dx  dr 

Equation  of  Radial  Motion 


d 

dt 


(mj  + 


-  (ouv)  + 
dx 


(1A)  ~  (ov2r) 

dr 


d_ 

dx 


(Txr)  +  (1/r) 


-  (ro,r)  -  aH 
dr 


(11) 
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Energy  Equation 


.  +  -  (one)  +  (1  /r)  '  (over)  =  — 

dt  dx  t>r  "  3x 


de  /  /.i  e 

yTx^ +  pO  “  +  "'■> 


+  (I/O 


d 

."Sr 


<lt’  .  u 
dr  K  Pr 


+  n  )  ~  r  (UT*r  +  V(Jrr) 

Pr, 


(12) 


These  equations  are  written  in  vector  form  as  (4:17): 


dU  dF 

- 4-  - 

dt  dx 


+  (1/r) 


dG 

dr 


=  (1/r)  H 


(13) 


where : 


0 


Q 

QU 

QV 

.  Q e  . 


F 


QU 

oir  - 
quv  -  r(, 

de  .  n  f  .  . 

Que  -  y  (—  +  — )  -  (mcj„  +  vr„) 
dx  Pr  Pr, 


Qvr 

quit  -  lur 
gv2r  -  onr 

de  .  u  e  . 

over  -  ry  - +  — )  -  r  (uzxr  +  vtjJ 
<V  Pr  Pr, 
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H 


'  0 
0 

-oH 

0 


The  transformed  vector  equation  for  general  curvilinear 
coordinates  is : 


dU  dF  (IF  «■)/; 

V  +  at  iii  +  7,  ,-ix  +  (l/r) 


dG  dG  drj 

d£  dr  drj  dr 


( 1  A)  H 


(  14) 


Applying  the  explicit  MacCormack  scheme  to  the  vector 
equation  results  in  the  following  algorithm  (4:27): 


{/(£  .  tj  .  /  +  A/)  =  L(A/)  A(C  .  7  .  r) 


(15) 


where  L(At)  is  the  two-dimensional  numerical  operator 
representing  MacCormack 's  algorithm  acting  on  the 
transformed  conservation  equations.  Through  use  of  time¬ 
splitting,  this  two-dimensional  operator  is  separated  into 
two  one-dimensional  sweep  operators  in  the  £  -  and  7  - 
directions.  The  operator  (At )  represents  the  solution  of 
the  equation: 


dU 

- + 

dt 


dP 

<l£  dx 


+  <>/'>§# 

di,  dr 


=  0 


(16) 


in  the  £  -  direction  by  a  time  increment  of  At  seconds.  In 
a  like  manner,  the  operator  L^(At)  represents  the  solution 
of  : 


dU  dF  drj 

dt  drj  dx 


+  d/r) 


d G  drj 
drj  dr 


=  dA)  H 


(17) 


in  the 


7  - 


direction 


by  a  time 


increment  of  A t  seconds. 
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The  dependent  variable  vector  U(  £  ,  rj  ,  t)  can  then  be 
advanced  in  time  as: 


U(Z  .>/./  +  At)  =  [  L',/2(A t/M)  Lv(At)  \"/2{At/M)\  U(Z  ,  >/  .  /)  (18) 

with  At  =  At*  if  A tK  <  A/, 
or  as : 


U(C  .  rj  .  /  +  At)  =  [  Vi/2(At/N)  L^Ar)  L?/2(A //AO]  .  V  •  0  (19) 


with  A/  =  A/,  if  Ar,  <  A tK 

where  M  and  N  are  the  smallest  even  integers  of  the 
quotients  (At?/Atc)  and  (At?  /At^  )  respectively.  The 
timesteps  At^  and  At^  are  the  maximum  allowable  timesteps  in 
the  C  -  and  ^  -  directions  as  determined  by  stability 
requirements . 

The  finite  difference  form  of  these  sweep  operators 
consist  of  a  two  step  predictor  -  corrector  procedure  which 
uses  one-sided  differencing  in  the  direction  of  sweep  and 
central  differencing  in  the  direction  perpendicular  to  the 
direction  of  sweep.  When  complete,  this  method  is 
equivalent  to  a  second  order  central  differencing  scheme  in 
two  dimensions. 
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The  L^(4t)  sweep  operator  represents  the  following 
numerical  procedure  (4:29): 

PREDICTOR 


(Fk.jT  -  (Fk-ijT\  (fxh.j 


-  (yr)~[ruj  (Gu.jT  ~  rk.,  (Gk.hjT]  (fyk.j 

(20) 


CORRECTOR 


(0k.,r+]  =  Vi  { (Ou.iT  +  ( uk.jT 


+  '/; 


-  (h,r 


+  '/: 


■i  <£),, 


-  (c,.,.,r'4  - m.;  (0t,/r*'/‘i  Sy, } 

^  (21) 

Similarly,  the  Lij(At)  operator  represents  the  following 
predictor  -  corrector  steps  (4:29,31): 

PREDICTOR 


(0u..,T  +  '/:  =  (0k.,T  -  ^  [(4,r  -  (Fu.j-iT\ 


-  (>A)£  K>  (^.,r  -  rk,  (Gk.hxT\  (Jk/  +  (1A)  A/  (A*./ 


(22) 
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CORRECTOR 


-  ,/;  <  ^'r  *  ‘ft./**  -  -  (\,r'i  A. 


A, 


av 


A/ 


(1A)^ir*.y  «^v,+  ir,;  -  A.;  iG*.,r+,i  (-)*.,  +  oa)  a/  (>/*./  } 


(23) 


where  n  and  h-*1  designate  values  at  time  t  and  t  +  At 
respectively  and  n+1/2  indicates  a  temporary,  predicted 
value.  The  values  ^ ,  —■ ,  ,  and  f2-  depend  only  on  the 

?x  3  r  3*  »r 

grid  geometry  and  are  calculated  only  once  for  each  grid 
point . 

Variants  of  this  program  for  two-dimensional  internal 
flow  and  three-dimensional  external  flow  have  been  used  in 
studies  by  Olson,  McGowan,  and  MacCormack  (10);  Shang, 
Hankey  and  Smith  (13);  and  Shang  and  MacCormack  (14).  The 
internal  axisymmetric  version  of  the  code  was  developed  for 
in-office  use  by  WRDC/FIMM. 


Incorporation  of  Nozzle  Structure 

In  order  to  model  the  ejector  geometry  desired,  it  was 
necessary  to  identify  boundary  conditions  which  accounted 
for  the  nozzle  structure  and  the  flow  injection  at  the 
nozzle  exit  in  addition  to  the  usual  wall  and  centerline 
conditions.  Once  the  desired  location  of  the  nozzle  exit 
was  selected,  the  grid  was  refined  in  that  area  to  better 
handle  the  large  gradients  which  exist  at  the  nozzle  exit. 
Then,  the  grid  points  which  best  approximated  the  corners  of 
the  desired  nozzle  exit  were  identified  (Figure  10). 
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Fig.  10.  Identification  of  Nozzle  Exit  Corners 


The  outside  walls  of  the  nozzle  structure  were  defined 
as  being  three  grid  lines  outside  the  nozzle  exit.  The 
outside  walls  and  ends  of  the  nozzle  structure  were  defined 
as  SOLID  WALLS  for  the  boundary  conditions  (Figure  11). 

The  grid  points  in  the  gap  were  defined  as  the  NOZZLE 
EXIT  for  the  boundary  conditions.  The  area  enclosed  within 
the  nozzle  structure  did  not  effect  the  computation  of  the 
flowfield  and  was  set  to  zero  velocity.  Figure  12  shows  a 
diagram  of  the  grid  boundaries  with  the  types  of  boundary 
conditions  applied  identified. 
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Fig.  11.  Identification  of  Nozzle  Walls 


Boundary..  Conditions 

Boundary  conditions  were  imposed  to  accommodate  the 
geometrical  restrictions  on  the  flowfield  by  defining  the 
solid  surfaces,  to  incorporate  primary  nozzle  flow  based  on 
a  constant  plenum  pressure,  and  to  provide  the  necessary 
upstream  and  downstream  conditions  to  achieve  closure. 
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Fig.  12.  Application  of  Boundary  Conditions 


Solid  Walls.  For  the  wall  of  the  ejector  and  the  inner 
and  outer  walls  of  the  primary  nozzle,  applying  the 
conditions  of  no  penetration  and  no  slip  provides: 

(Q“)w  =0  (24) 

(Qv\*  =0  (25) 

3  p 

Allowing  no  normal  pressure  gradient  (  —  =  0  )  leads  to: 

an 

(QL’)w  =  Q\  [  e\-Vt  (u\  +  vf)  ]  (26) 

Applying  a  constant  wall  temperature  T»: 

(Q)~  =  (Qe)w  /  (cvT„)  (27) 

where  the  subscript  w  indicates  the  value  at  the  wall  and 
the  subscript  1  indicates  a  value  taken  one  grid  point  away 
from  the  wall. 
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Center  I  ine .  Since  the  centerline  is  the  axis  of 
symmetry,  there  can  be  no  cross  velocity.  This  results  in 
the  equation: 


(<w)(7_  =  0  (28) 

Also,  there  can  be  no  radial  gradients  of  u,  p,  or  T  (  r  : 

3  e  3T  ^  aT( 

=  —  =  0  ),  resulting  in  the  following  equations: 

L>cl  =  Q  l  (29) 

=  (QU)  1  (30) 

(Qe)cL  =  ((?<?) i  -  '/2  {Qvfx  /  Qcl  (31) 

where  the  subscript  CL  indicates  the  value  at  the  centerline 
and  the  subscript  1  indicates  a  value  taken  one  grid  point 
away  from  the  centerline. 

Downstream .  Enforcing  conservation  of  mass  in  the  £  - 

direction  (  37^  +  1/r  77—-  =  0  )  by  allowing  -  0  and 

>-&■  -  0  gives:  {  K 

<>  s 


(Q“)k 

=  (Qu)k-\ 

(32) 

(Qv)k 

=  (<?w) K-l  /  rK 

(33) 

where  the  subscript  K  indicates  the  value  at  the  downstream 
boundary  and  the  subscript  K-l  indicates  a  value  taken  one 
grid  point  upstream.  Since  the  flow  in  the  ejector  modeled 
is  exclusively  subsonic,  the  exit  static  pressure  must  equal 
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the  ambient  pressure  Pa.  Applying  this  constraint  and 
allowing  temperature  to  propagate  downstream  leads  to: 

(l>)k  -  Pa  /  RTk- i  (34) 

(<*)*  =  Ok  [  {cJk- i)  +  Vt  (u\  +  \%)  |  (35) 


Upstream .  Enforcing  conservation  of  mass  in  the  C  - 
direction  as  was  done  on  the  downstream  boundary  leads  to: 


({?«)  i 

=  (e«)2 

(36) 

(nv), 

=  {Qvr)2  /  rx 

(37) 

where  the  subscript  1  indicates  the  value  at  the  upstream 
boundary  anu  *  .,e  subscript  2  indicates  a  value  taken  one 
grid  poir.  .uwnstream.  Additionally,  assuming  incom- 
pressitie,  isentropic  flow  from  ambient  stagnation 
conditions  results  in  the  equations: 

(<?)i  =  Pa  /  RTa  (  38  ) 

(C><? )l  =  Ch  [  {cJu)  +  Vi  («i  +  vf)  (y-  1  )/y  ]  (39) 

where  Pa  and  T*  are  the  ambient  stagnation  pressure  and 
temperature  respectively. 

Nozzle  Exit.  The  nozzle  exit  was  treated  the  same  as 
the  upstream  boundary  except  that  the  nozzle  total  pressure 
(Po  no*)  was  based  on  the  defined  pressure  ratio  (Prat): 

Pt)noz  =  Pral{Pa)  (40) 

and  the  direction  of  the  flow  was  constrained  to  the  desired 
injection  angle  (<Xi). 
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Turbulence  model 


The  turbulence  model  used  in  the  original  program  was  a 
modified  Baldwin-Lomax  model  based  on  the  upper  surface. 

This  turbulence  model  uses  a  two  region  calculation  of  the 
eddy  viscosity  (14:4).  The  inner  region  calculation  f 1  is 
used  up  to  the  distance  y*  from  the  wall  after  which  the 
outer  region  formulation  to  is  used.  The  distance  y*  is  the 
defined  as  the  smallest  distance  from  the  wall  at  which  £1 
and  £o  are  equal. 

The  inner  formulation  is  given  by  (14:4): 

f;  =  e  (0-4  L  D)2  I"  I  (41) 

where  the  vorticity  of  the  flow,  60,  the  Van  Driest  damping 
factor,  D,  and  the  scaling  length,  L,  are  given  by: 


=  Vi  V  x  U 

(42) 

I)  =  \  -  exp[  {  -  qw  \(ow  |  /  ftw  )Vl  L  /  26  | 

(43) 

L  =  [  (jt  -  xw)2  +  (r  -  rw)2  )'/; 

(44) 

and  the  subscript  w  represents  the  value  at  the 

outer  formulation  is  given  by  (14:4): 

wall.  The 

=  0.0336  q  Fwake  /  [  1  +  5.5  (0.3  L/LmMf  ] 

(45) 

where  the  wake  function,  Fwake,  is  the  minimum 

following  at  any  point  in  space: 

value  of  the 

F wake  ~  Lm3X  Fmax 

or 

Fwake  ~  0.25  Lmax  Kmax  /  ^max 

(46) 
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and  Ln&x  is  the  value  of  the  length  scale,  L,  where  F: 

F  =  L  D  \ui  |  (47) 

reaches  its  maximum  value,  Foax,  within  the  turbulent  shear 
layer . 

To  simplify  this  investigation,  no  attempt  was  made  to 
identify  better  turbulence  models  for  this  type  of  wall  jet 
flow.  The  turbulence  model  already  incorporated  in  JSIAXK 
was  not  changed  but  was  based  on  the  nearest  wall  to  include 
nozzle  wall  effects. 


Er..Q£.£dur.s. 

Once  the  desired  location  of  the  nozzle  exit  was 
identified,  the  vertical  and  horizontal  grid  spacing  tables, 
EJECVRT  and  EJECHRZ  (Appendices  A-3.A-4),  were  modified  to 
refine  the  grid  in  that  vicinity.  EJECGRD  was  then  run  to 
produce  the  grid  and  establish  the  initial  flowfield. 

Setting  the  flow  injection  angle  (<Xi )  to  approximate  the 
flow  direction  established  by  the  simplified  initial 
flowfield  (to  minimize  extreme  gradients)  the  flow  solver 
JSIAXK  could  be  started. 

Initially,  the  timestep  per  iteration  was  kept  small 
until  the  larger  gradients  could  be  smoothed  out,  then  was 
increased  to  provide  faster  convergence  to  the  steady  state. 
Once  the  large  gradients  of  the  initial  flowfield  were 
smoothed  out,  it  was  also  possible  to  change  di  until  the 
desired  value  was  set.  The  program  was  then  run  until 
steady  state  was  reached.  Steady  state  was  defined  for  the 
purpose  of  this  investigation  as  a  change  of  less  than  0.001 
ft/sec  in  total  velocity  at  each  of  four  critical  locations 
in  3000  iterations.  This  process  took  approximately  90,000 
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iterations  to  complete  at  0.151  CPU  seconds  per  iteration  on 
the  CRAY  XMP  computer.  This  solution  was  then  stored  on  the 
Central  Datafile  System.  A  new  was  then  set  to  the  next 
angle  of  interest.  Convergence  on  a  new  steady  state 
solution  then  took  approximately  3600  iterations. 

This  process  was  continued  until  solutions  were 
achieved  for  all  desired  injection  angles  for  the 
conf igurat ion  under  investigation.  For  a  new  nozzle 
location,  the  process  was  restarted  from  the  beginning. 
Configurations  la,  lb,  and  lc  were  able  to  be  run 
consecutively  because  the  changes  to  the  flowfield  were 
minimal  in  switching  from  one  to  the  other. 

The  post-processing  program  AUGMENT  was  then  run  for 
each  flowfield  solution.  This  program  integrates  the  x  - 
momentum  across  the  exit  and  compares  it  to  the  isentropic 
thrust  of  the  nozzle  alone  to  calculate  the  thrust 
augmentation  per  the  control  volume  analysis  described  in 
section  I.  A  sample  comparison  of  mass  flows  in  and  out  was 
also  accomplished  to  ensure  that  conservation  of  mass  was 
being  satisfied. 
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Ill 


Results  and  Discussion 


Case.s  Investigated 

The  cases  investigated  consist  of  four  nozzle  location 
configurations  fc  one  ejector  geometry.  The  geometry  of 
the  ejector  was  .elected  to  be  similar  to  one  investigated 
by  Unnever  to  allow  for  direct  comparison.  The  dimensions 
of  the  ejector  studied  by  Unnever  (Figure  7)  were  (17:60): 


Li  = 

3 . 5 

in 

Yi  =  3° 

1.2  = 

3 . 5 

in 

N 

1! 

CD 

0 

resulting  in  a  diffu 

ser 

area  ratio 

(y9) 

of  1.71. 

The  dimensions 

of 

the  annular 

nozz 

le  of  Unn 

experiment  and  those  of  the  configurations  investigated  are 
tabulated  in  Table  I.  The  dimensions  of  Unnever 's  nozzle 
indicate  the  it  was  0.55  inches  from  the  wall.  Unnever 's 
testing  was  carried  out  with  a  pressure  ratio  of  1.14  which 
was  replicated  for  all  configurations  investigated. 

Configuration  la.  The  location  and  geometry  of  the 
nozzle  for  configuration  1  was  chosen  to  approximate  as 
closely  as  possible  that  of  the  annular  nozzle  tested  by 
Unnever . 

Configurations  0.  2.  3.  The  locations  and  geometries 
of  the  nozzles  for  configurations  0,  2,  and  3  were  chosen 
such  that  the  distance  from  the  wall  and  the  nozzle  area 
would  be  approximately  the  same  as  configuration  la.  This 
allowed  investigation  of  the  thrust  augmentation  effects  of 
nozzle  position  angle,  8,  only  (Figure  13). 
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TABLE  I 

Nozzle  Geometries 


Li 

Radius 

Gap 

Aj 

0 

a 

i 

;  Unnever 

1 

2.25 

3.00 

0.065 

1.226 

61 . 9 

0.081 

1 

!  Config  la 

2.25 

2 . 96 

0.072 

1.349 

61.9 

0.089 

|  Config  0 

2.40 

3.30 

0.064 

1.325 

70.0 

0.087 

Config  2 

1.96 

2.51 

0.081 

1.296 

50.0 

0.085 

;  Config  3 

1 . 64 

2 . 17 

0.094 

1.306 

40.0 

0.086 

Config  lb 

2 . 26 

2.96 

0.053 

0.989 

61.9 

0.065 

Config  lc 

2.25 

2.95 

0.095 

1.774 

61.9 

0. 117 

Configurations  lb.  lc .  The  locations  of  the  nozzles 
for  configurations  lb  and  lc  were  chosen  to  be  as  close  to 
that  of  configuration  la  as  possible.  The  only  change  made 
was  to  the  size  of  the  gap  and  therefore  the  nozzle  area. 
This  was  done  to  investigate  nozzle  area  ratio  effects  on 
thrust  augmentation. 


Primary  .flbjaa.tiY.ajs 

Comparison  with  Experimental  Data.  Configuration  la 
was  designed  to  model  as  closely  as  possible  the  annular 
nozzle  configuration  used  in  Unnever's  experimental 
investigation.  A  qualitative  comparison  of  the  exit 
velocity  profile  of  configuration  la  with  an  experimentally 
measured  profile  from  Reznick  (Figs.  14  and  15)  shows  that 
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Fig.  14.  Exit  Velocity  Profile 


Numerical 


Fig.  15.  Exit  Velocity  Profile  -  Experimental  (11:48) 
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the  shape  of  the  profile  is  closely  replicated  by  the  model. 
This  is  an  indication  that  the  flow  structure  is  indeed 
similar  to  that  in  the  experimental  ejector.  However,  since 
the  experimental  data  was  produced  using  the  hook  nozzles 
and  is  therefore  only  similar  to  the  ejector  modeled,  the 
comparison  cannot  be  exact. 

Similarly,  comparison  of  thrust  augmentation  vs. 
injection  angle  curves  for  the  numerical  model  (Figure  16) 
with  similar  curves  from  Unnever  (Figure  17)  (17:46)  show 
that  the  expected  trends  are  indeed  modeled.  Again,  since 
the  experimental  curves  are  from  hook  nozzle  cases  (by 
necessity,  since  the  injection  angle  of  the  annular  nozzle 
could  not  be  varied)  the  comparison  can  only  be  qualitative, 
not  quantitative. 


Fig.  16.  Thrust  Augmentation  vs.  Injection  Angle 

Numerical 
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Fig.  17.  Thrust  Augmentation  vs.  Injection  Angle 

Experimental 


However,  a  quantitative  comparison  of  the  thrust 
augmentation  from  configuration  la  with  Unnever's  annular 
nozzle  configuration  (Table  II)  (17:60)  shows  that  the 
numerical  model  did  not  achieve  the  thrust  augmentation 
measured  experimentally. 


Table  II 

Thrust  Augmentation  Comparison 


Ej  ector 

0 

Numerical  Model 

Experimental 

1.40 

1.48 
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Possible  explanations  for  this  difference  include: 


1.  Turbulence  modeling  not  optimized  for  a  wall  jet 
flow  resulting  in  poorer  mixing  in  the  numerical  model. 

2.  Three  dimensional  mixing  or  swirl  in  the 
experimental  ejector  which  is  not  accounted  for  in  the 
numerical  model  due  to  the  assumption  of  axial  symmetry. 

3.  Modeling  the  jet  exiting  the  nozzle  as  perfectly 
parallel  flow  in  the  desired  direction  when  there  is  likely 
some  divergence  in  the  actual  jet  experimentally. 

Each  of  these  lead  to  less  complete  mixing  in  the 
numerical  model  and  would  thus  account  for  the  lower  thrust 
augmentation  calculated.  Additionally,  other  sources  of 
error  which  would  account  for  differences  in  experimental 
versus  numerical  results  include: 

4.  Insufficient  grid  resolution. 

5.  Experimental  error. 

Using  0  -  1.00,  or  no  augmentation,  as  a  baseline,  this 
results  in  a  16.7  percent  error.  This  error  indicates  that 
the  model,  as  it  is,  would  not  be  a  good  predictor  of  the 
actual  thrust  augmentation  which  could  be  expected  from  an 
experimental  setup.  However,  since  trends  determined  with 
this  model  appear  to  be  consistent  with  experimental 
results,  analysis  of  different  parameters  should  allow 
identification  of  values  which  maximize  thrust  augmentation 
which  could  be  used  to  guide  hardware  development. 

Data  Visualization.  One  of  the  benefits  of  a  numerical 
simulation  is  that  the  solution  results  in  complete 
knowledge  of  the  computed  flowfield.  This  knowledge  can  be 
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used  to  provide  insight  into  the  behavior  of  the  flow  within 
the  ejector  which  would  be  nearly  impossible  to  determine 
experimentally.  Plotting  velocity  contours  (Figure  18) 
clearly  shows  the  primary  characteristics  of  the  flowfield. 


V/Uref  Uref  =  1U-7  ft/sec 
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Fig.  18.  Total  Velocity  Contours 

First,  (noting  that  the  flow  is  from  left  to  right)  one  can 
see  how  the  jet  hugs  the  wall  and  negotiates  the  curvature 
in  the  inlet  due  to  the  coanda  effect.  Second,  one  can  see 
how  the  jet  spreads  out  as  it  entrains  more  of  the 
surrounding  flow  along  with  it.  The  increasing  velocity  of 
the  entire  flowfield  as  the  cross-sectional  area  of  the  duct 
decreases  can  also  be  seen.  Lastly,  a  relatively  high 
velocity  channel  flow  can  be  seen  between  the  nozzle 
structure  (and  the  jet)  and  the  wall. 

A  streamline  diagram  of  configuration  la  at  its  optimum 
injection  angle  (Figure  19)  shows  the  smooth  flow  which 
characterizes  the  best  performing  injection  angles. 
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Fig.  19.  Streamlines 


One  can  also  sea,  however,  that  a  stagnation  point 
exists  on  the  nozzle  structure  upstream.  Although  the 
annular  nozzle  structure,  as  modeled,  is  thinner  than  that 
of  the  experimental  apparatus,  this  diagram  shows  that 
significant  flow  blockage  by  the  nozzle  exists  even  in  the 
numerical  model  just  as  was  suspected  in  the  experimental 
work  by  both  Reznick  and  Unnever.  A  velocity  vector  plot  of 
this  region  (Figure  20)  clarifies  the  picture. 

Finally,  a  sequence  of  plots  of  velocity  profiles 
(Figures  21-23)  shows  how  the  jet  mixes  with  the  secondary 
flow  as  it  progresses  through  the  ejector.  The  first  plot 
(Figure  21)  shows  the  velocity  profile  at  three  grid 
stations  ranging  from  the  exit  plane  of  the  nozzle  (K  =  32) 
to  a  cross-section  within  the  constant  area  mixing  chamber, 
or  throat,  of  the  ejector  (K  =  75).  These  plots  show 
clearly  how  the  jet  begins  to  mix  with  the  secondary  flow 
and  how  the  entire  flow  accelerates  as  the  inlet  converges. 
It  may  further  be  determined  from  this  plot  is  that  the  grid 
resolution  at  the  upper  surface  appears  to  be  barely 
sufficient  to  capture  the  boundary  layer.  However,  it  does 
not  appear  to  be  sufficient  to  capture  boundary  layer 
effects  on  the  nozzle  walls. 
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The  second  plot  (Figure  22)  shows  how  the  jet  continues 
to  mix  with  the  secondary  flow  as  it  progresses  from  the 
beginning  (K  =  75)  to  the  end  (K  =  100)  of  the  mixing 
chamber . 


Fig.  22.  Mixing  Chamber  Velocity  Profiles 


The  last  plot  (Figure  23)  portrays  the  changing 
profiles  as  the  flow  proceeds  through  the  diffuser  from  the 
end  of  the  mixing  chamber  (K  =  100)  to  the  exit  (K  =  130). 
The  deceleration  of  the  flow  in  the  diffuser  is  evident. 
Also,  the  continued  spreading  of  the  jet  as  it  mixes  with 
the  secondary  flow  is  clearly  visible.  These  basic  flow 
features  were  consistent  for  all  of  the  configurations 
modeled . 
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Fig.  23.  Diffuser  Velocity  Profiles 


Flexibility .  This  code,  as  modified  for  this 
investigation,  is  capable  of  handling  variations  of 
injection  angle,  nozzle  position,  nozzle  gap  (and  therefore 
area)  and  pressure  ratio.  With  some  additional  effort, 
diffuser  and  inlet  geometries  could  be  altered  as  well. 

Secondary  Objectives 

Several  of  the  secondary  objectives  were  approached  in 
the  course  of  tnis  investigation  although  all  would  require 
further  work  to  provide  the  requisite  parametrics  to  use  for 
hardware  design.  The  parameters  studied  include  injector 
angle,  nozzle  position,  and  nozzle  to  throat  area  ratio. 
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In.iection  Angle.  The  effect  of  injection  angle  on  the 
flowfield  can  be  observed  by  comparing  the  optimum  injection 
angle  for  configuration  la  to  two  off-nominal  cases.  For 
configuration  la,  the  optimum  injection  angle  was  found  to 
be  «•  i  =  52  degrees.  This  is  compared  to  the  off-nominal 
cases  of  cti  =  64  degrees  and  =  36  degrees.  The  thrust 
augmentation  achieved  by  each  of  these  cases  is  tabulated 
in  Table  III. 


Table  III 


Augmentation  with  Changing  Injection 


Angle  (<*i) 


*1 

0  | 

36 

1.264  j 

52 

1.398  | 

64 

1.308 

First,  looking  at  velocity  contours  for  the  =  64  degree 
case  (Figure  24)  one  can  see  that  the  jet  flow  does  not 
follow  the  wall  closely  but  instead  spreads  out  into  the 
main  mixing  chamber.  The  improved  mixing  accomplished  by 
the  increased  injection  angle  would  seem  to  be  desirable. 
But,  as  the  exit  velocity  profiles  indicate  (Figure  25),  the 
increased  centerline  flow  is  accomplished  at  the  expense  of 
a  large  drop  in  flow  velocity  near  the  wall.  For  a  circular 
ejector,  the  region  near  the  wall  represents  a  much  greater 
area  in  which  momentum  was  sacrificed  than  the  area  near  the 
centerline  where  momentum  was  gained.  Also,  it  can  be  seen 
that  the  flow  has  begun  to  separate  from  the  wall  resulting 
in  reversed  flow. 
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Fig.  24.  Velocity  Contours  :  <*- 1 


64  degrees 


The  second  off-nominal  case  (  a ±  =  36  degrees)  is  one 
in  which  the  jet  is  directed  too  near  the  wall.  A  first 
look  at  the  velocity  contours  (Figure  26)  shows  the  jet 
closely  hugging  the  wall  which  by  the  above  argument  would 
seem  to  be  desirable.  However,  a  closer  look  at  the  area 
between  the  jet  and  the  wall  shows  that,  contrary  to  the 
optimum  case  (Figure  18),  the  flow  through  this  channel  is 
of  relatively  low  velocity.  Also,  the  high  speed  jet  flow 
along  the  wall  loses  more  momentum  due  to  frictional 
effects . 


The  effect  of  the  blockage  of  this  high  energy  flow 
source  and  frictional  losses  are  reflected  in  the  exit 
profile  (Figure  27).  Although  the  jet  is  more  concentrated 
for  the  i  =  36  degree  case  it  is  no  more  powerful  than  the 
optimal,  a.  i  =  52  degree,  case.  The  reduced  mixing  also 
results  in  a  decrease  in  momentum  across  the  majority  of  the 


center  section  of  the  ejector.  Apparently,  the  optimum  case 
is  one  which  balances  these  considerations  by  achieving 
maximum  mixing  without  allowing  the  jet  to  separate  from  the 
wall. 


Fig.  27.  Exit  Velocity  Profiles  :  i  =  52  and  36  degrees 


Figure  28  defines  a  pair  of  angles  which  the  injection 
angle  of  the  nozzle  can  be  referenced  to.  These  angles  are 
those  representing  injection  parallel  and  tangent  to  the 
wall.  Figure  29  shows  the  relationship  of  the  injection 
angle  for  maximum  thrust  augmentation  to  those  representing 
injection  parallel  to  the  wall  and  tangent  to  the  wall  for 
each  of  the  four  configurations  investigated.  In  all  cases, 
the  optimum  injection  angle  is  approximately  a  quarter  of 
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the  way  between  parallel  and  tangent  to  the  wall.  This 
conclusion  must  be  caveated  by  noting  that,  due  to  the 
circular  inlet  lip  and  the  fact  that  all  configurations  were 
limited  to  be  the  same  distance  from  the  wall,  the  relative 
geometry  between  the  nozzle  and  the  wall  is  the  same  for  all 
configurations.  To  fully  investigate  this  parameter, 
additional  data  would  need  to  be  collected  at  different 
distances  from  the  wall  and  (if  desired)  for  different  inlet 
geometries . 


Fig.  28.  Reference  Injection  Angles 


Nozzle  Position.  In  investigating  the  influence  of 
nozzle  position  on  thrust  augmentation,  it  was  necessary  to 
isolate  the  effect  of  nozzle  movement  by  constraining  the 
conf igurat  ions  chosen.  In  order  to  eliminate  changing  wall 
effects,  all  configurations  were  set  to  have  the  nozzle  the 
same  distance  from  the  wall  as  configuration  la  (0.55  in.). 
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Fig.  29.  Optimum  Injection  Angle  with  Relation 
to  Reference  Angles 


The  nozzle  gap  size  was  adjusted  to  achieve  similar  nozzle 
areas  for  all  configurations  to  eliminate  area  ratio 
effects.  This  reduces  the  effect  of  location  changes  to 
only  the  effect  of  nozzle  position  angle,  0.  Plotting  the 
maximum  augmentation  achieved  at  each  location  vs.  the 
posit. ion  angle  (Figure  30)  a  definite  trend  can  be  detected. 
Higher  augmentation  ratios  are  achieved  with  decreasing 
position  angle,  or  in  other  words,  placing  the  nozzle 
further  into  the  inlet. 

The  effect  of  nozzle  location  can  best  be  seen  by 
comparing  the  two  extremes  studied,  configuration  0  and 
configuration  3.  By  looking  at  velocity  profiles  of  the  jet 
itself  (Figures  31  and  32)  we  can  see  that  configuration  3 
has  a  higher  peak  jet  velocity  than  configuration  0. 
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Fig.  30.  Augmentation  vs.  Position  Angle 


Fig.  31.  Nozzle  Velocity  Profile  :  Configuration  3 
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Fig.  32.  No22le  Velocity  Profile  :  Configuration  0 


This  is  due  to  the  fact  that  the  plenum  pressure  of 
both  configurations  is  the  same: 

Po  noz  -  Prat  (Pa)  (47) 

but,  as  can  be  seen  in  Figures  33  and  34,  the  exit  of 
configuration  3's  nozzle  is  in  the  lower  pressure  region 
nearer  the  throat  of  the  ejector.  The  pressure  near  the 
exit  of  configuration  3's  nozzle  (Figure  33)  is  about  0.990 
Pa  on  the  side  nearer  the  wall  and  about  0.994  Pa  on  the 
opposite  side.  This  is  compared  to  configuration  Os  nozzle 
(Figure  34)  which  is  located  where  the  pressures  are  approx¬ 
imately  0.996  Pa  and  Pa  respectively.  Since  the  exit 
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pressure  is  lower  for  configuration  3,  the  exit  velocity  is 
necessarily  higher.  Also,  there  is  less  momentum  loss  due 
to  friction.  This  model  does  not,  however,  account  for 
losses  within  the  nozzle  which  would  reduce  the  effec¬ 
tiveness  of  the  longer  nozzle  structure  of  configuration  3. 
One  would  assume  that  there  would  be  some  point  at  which  the 
advantages  of  moving  the  nozzle  toward  the  throat  would  be 
outweighed  by  the  reduced  mixing  length  allowed.  Further 
analysis  would  be  necessary  to  define  the  optimum  location. 
One  can  also  see  in  these  plots  how  a  pressure  difference  is 
maintained  across  the  jet  it  which  forces  the  jet  to  turn. 
This  is  similar  to  the  turning  of  the  flow  from  a  jet  flap. 

Nozzle  Area  Ratio.  Configurations  lb  and  lc  were 
designed  to  keep  all  other  factors  the  same  as  configuration 
la  except  the  nozzle  gap,  and  therefore  the  nozzle  area,  in 
order  to  investigate  the  effect  of  nozzle  to  throat  area 
ratio  changes.  Configuration  lb  had  a  reduced  nozzle  area 
and  configuration  lc  had  an  increased  nozzle  area  relative 
to  configuration  la.  The  area  ratios  and  achieved  thrust 
augmentation  were  as  in  Table  IV.  The  results  are  plotted 
in  Figure  35. 


Table  IV 

Augmentation  with  Changing  Nozzle  Area  Ratio  (cx ) 


Conf ig . 

l  .  ■  .  - 

— 

Nozzle  Area  Ratio 

0 

'  lb 

0 . 0650 

1.407 

la 

0.0887 

1 . 398 

lc 

_ 

0. 1167 

1 . 353 
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Fig.  35.  Augmentation  vs.  Nozzle  Area  Ratio  (a) 


As  can  be  seen  from  Figure  35,  further  reduction  in 
area  ratio  would  appear  to  be  fruitless,  contrary  to  the 
theoretical  analysis  presented  earlier  which  predicted  ever 
increasing  augmentation  with  decreased  area  ratio.  But,  as 
Figure  36  shows,  these  results  are  in  line  with  experimental 
results  from  Fought's  MS  Thesis  as  published  in  McCormick. 
These  results  suggest  that  the  area  ratio  used  by  Reznick 
and  Unnever  (equivalent  to  configuration  la)  is  near  the 
optimum  for  this  size  ejector  (Diameter  =  112  mm). 
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Fig.  36.  Augmentation  vs.  Nozzle  Area  Ratio  ( l/oO 
Numerical  and  Experimental  (9:285) 
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IV .  Cane lus ions 


All  primary  objectives  of  this  thesis  have  been  been 
met.  The  numerical  model  has  been  shown  to  provide  a 
reasonable  simulation  of  the  flow  characteristics  through 
the  selected  ejector  configuration.  Although  the  simulation 
is  not  sufficient  to  permit  accurate  prediction  of  actual 
thrust  augmentation  values  for  experimental  hardware,  all 
trends  appear  to  be  adequately  modeled.  Taking  advantage  of 
the  complete  knowledge  of  the  flowfield  provided  by  a 
numerical  solution,  the  dynamics  of  the  computed  internal 
flow  were  able  to  be  analyzed.  Data  visualization 
techniques  were  able  to  give  us  insight  into  the  internal 
flow  of  the  ejector  which  could  not  be  measured 
experimentally.  Sufficient  flexibility  has  been 
incorporated  in  the  code  to  handle  variations  of  parameters 
such  as  injection  angle,  nozzle  position,  nozzle  to  throat 
area  ratio,  pressure  ratio,  as  well  as  inlet  and  diffuser 
geometries . 

In  addition,  the  secondary  objective  of  exploring  the 
effect  of  these  parameters  on  thrust  augmentation  was 
partially  investigated  allowing  some  initial  conclusions  to 
be  drawn  which  could  affect  future  experimental  hardware 
design . 

1)  The  area  ratio  used  by  Reznick  and  Unnever  is  near 
the  optimum  and  should  be  sufficient  for  future  experiments. 
However,  if  it  is  within  manufacturing  capability,  some 
reduction  of  nozzle  area  would  be  beneficial. 

2)  Future  experiments  should  concentrate  on 
configurations  which  inject  the  jet  at  a  location  closer  to 
the  entrance  of  the  constant  area  mixing  chamber  (smaller  9) 
if  possible. 
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that  an  annular  nozzle 


3)  Initial  indications  are 
design  should  attempt  to  provide  an  injection 
approximately  1/4  of  the  way  between  parallel  to 
and  tangent  to  the  wall  for  the  selected  nozzle 


angle 
the  wall 
location. 
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V .  Suggestions  and  Recommendations 


Any  continued  research  in  this  area  should  concentrate 
on  expanding  the  parametric  analysis  in  the  areas  of  nozzle 
location  as  well  as  inlet  and  diffuser  geometries.  Nozzle 
locations  further  into  the  inlet  and  variations  in  distance 
from  the  wall  should  be  investigated.  Larger  diffuser 
ratios  should  also  provide  increased  augmentation  capability 
and  different  inlet  shapes  could  optimize  the  incoming  flow. 
Further  analysis  of  injection  angle  effects  for  the  new 
geometries  may  also  be  necessary. 

Another  primary  area  of  research  which  would  be 
beneficial  is  application  of  improved  turbulence  modeling, 
better  grid  refinement,  and  second-order  boundary  conditions 
to  improve  the  numerical  model  and  try  to  bring  the 
predicted  thrust  augmentation  levels  in  line  with 
experimental  results.  It  may  even  be  desirable  to  change  to 
a  different  solution  algorithm  to  reduce  computation  time. 

Lastly,  using  the  parametric  results  presented  here  to 
design  new  experimental  hardware  should  result  in  improved 
performance  and  would  serve  to  better  validate  this  model. 
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FORW«RO  1  BACKWARD  DIFFERENCING  FOR  PREDICTOR  t  CORRECTOR  r  GENERATE  THE  HEAT  FLUX  TERMS 
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FORWARD  &  BACKWARD  DIFFERENCING  FOR  PREDICTOR  i  CORRECTOR 
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